Josephson effects in one-dimensional supersolids 
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We demonstrate that superflow past an obstacle is possible in a solid phase in the one-dimensional 
Gross-Pitaevskii equation with a finite-range two-body interaction. The phenomenon we find is 
analogous to the DC Josephson effect in superconductors and we deduce the "Josephson relation 
" between the current and phase difi'erence of the condensates separated by the obstacle. We also 
discuss persistent current and nonclassical rotational inertia in annular container with a penetrable 
potential barrier. The phase diagram in the plane of the current and the interaction strength is given. 
Our result provides a simple theoretical example of supersolidity in the presence of an obstacle. 

PACS numbers: 67.80.-s 



I. INTRODUCTION 

Supersolid is the quantum solid with a superflow prop- 
erty. Although the possibility of supersolid has been dis- 
cussed theoretically about 40 years ago^"— , extensive and 
intensive studies have started since Kim and Chan have 
reported^"— the first evidence for the non-classical rota- 
tional inertiap^ in a solid phase of *He in the torsional 
oscillator experiments. Subsequently, experimental stud- 
ies have revealed various puzzling nature in supersolid 
■^He^^V'! enhancement of elasticity having the same tem- 
perature dependence as the resonance frequency shifti^, 
the sensitivity of resonance frequency shift on history of 
annealing processi^ and small concentration of "^He^^, 
superflow along the grain boundarieai^, ultraslow re- 
laxation dynamics in the dissipation and resonance fre- 
quency shift in the torsional oscillator experimental^. 
As consistent explanations of these experimental re- 
sults, several theories have been proposed such as defect- 
mediated superfluidity mechanismsi^iii, a non-superfluid 
glass picturei^ii^ and intrinsic supersolidity of defect-free 
crystaP^. 

Supersolid is also related to other systems; experi- 
ments on ^He film realized on graphite suggest a pos- 
sibility of two-dimensional supersolid^^. Bose-Einstein 
condensate(BEC) of ^^Cr might be another possible can- 
didate of the system with a supersolid phase, owing to 
the long-range dipole-dipole interactions22r^. The pos- 
sibility of supersolid of Rydberg atoms has been studied 
very recentl y^^'^^ . Thus, the properties of intrinsic su- 
persolids should be clarified in more detail from broader 
contexts. 

As a phenomenological model of supersolid, Pomeau 
and Rica2^ have investigated a modified version of the 
Gross-Pitaevskii (GP) modeP^''^°, which has a two-body 
repulsive interaction not of a contact-type but with a 
finite-range. The liquid has the excitation energy with 
roton minimum for sufficiently large repulsive interac- 
tion (or sufficiently high density of particles) . For larger 
repulsive interaction, the ground state acquires the po- 
sitional order detected by the Bragg peaks. It has been 



shown that the solid exhibits the non-classical rotational 
inerti a^^'^^" — and quantized vortices^S. under a uniform 
rotation. While the solid phase could be identified with a 
supersolid from these results, it was found that the dissi- 
pationless flow around an obstacle is not possible^* . It is 
thus important to see whether the supersolid can sustain 
the superflow in the presence of obstacles or not. This is 
the issue addressed in the present paper. 

As a simplest set-up, we consider the GP equation with 
a finite-range two-body interaction in one-dimensional 
systems. We show that the steady flow state exists 
in the solid phase even in the presence of an obstacle. 
The superflow in the presence of an obstacle in a one- 
dimensional system can be regarded as a DC Josephson 
current. That is, we examine superflow property of solid 
phase(, which we call supersolidity in this paper) through 
the DC Josephson effeclj^. The Josephson effects ob- 
served in ■*IIe and in BEG in cold atoms serve an evidence 
for superfluidity'^^-'^^. In addition to the Josephson effect, 
we also discuss two other important superflow proper- 
ties: the persistent-current effect and Hess-Fairbank ef- 
fect in supersolid state in a ring with a potential barrier. 
The former effect is a flow property of metastable states 
while the latter is a property of thermodynamically sta- 
ble state in which the whole or a part of supersolid is 
at rest against in a container rotating with a sufficiently 
small angular velocity. These two effects are considered 
as two fundamental properties of superfluids and cor- 
respond, respectively, to persistent electric current and 
Meissner effects in superconductors. 

This paper is organized as follows. In Sec. |TT1 we in- 
troduce the GP equation with a finite-range two-body 
interaction and explain the method of numerical calcu- 
lations. In Sec. mil we first present the phase diagram 
in the plane of current J and (dimensionless) two-body 
interaction strength g. We next present the results re- 
lated to the Josephson effects in the solid phase. Fur- 
ther we discuss persistent current and present results on 
nonclassical rotational inertia in Hess-Fairbank effect. In 
Sec. IIV| we discuss implications to systems in higher di- 
mensions and possible difference of superfluid properties 



between the supersolid phase and conventional superflu- 
ids. The conclusion is given in Sec. |Vl 



II. MODEL AND METHOD 

We start with a one-dimensional GP equation with a 
finite-range interaction^^ 

,d^(x,t) n^ d"^ ^, , ,^^, , 1 ,, , 

rL/2 
+ dyV{x-y)\^{y,t)\H{x,t)=0, (1) 

J-L/2 

where m is the atomic mass, ^(a;) is the condensate wave 
function. L is the system size. V{x) and U{x) represent, 
respectively, two-body repulsion and the potential bar- 
rier, /i denotes the chemical potential, which is deter- 
mined by the condition on the total number of particles 
being^l 



rL/2 
N= dx\'^{x)f 

J-L/2 



(2) 



In the presence of potential barrier U{x), there are two 
conservation laws^S,; one is equation of continuity of par- 
ticle density 
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with the current^ 
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and the other is equation of continuity of local energy, 
the latter of which we do not discuss in this paper. 

If the stationary solution {d'^/dt = 0) to ([Ij exists, 
this solution has the current J{x, t) being spatially con- 
stant (In the following, we denote by J a constant cur- 
rent ). The stationary solution with finite J represents 
a condensate with macroscopic dissipationless flow. The 
existence of such a stationary solution in the presence of 
potential harrier U{x) can be regarded as an evidence for 
superfluidity or supersolidity. In earlier works^^"— , sta- 
tionary solutions with finite J in the presence of a short- 
range barrier U{x) ex 5{x) were found when two-body 
interaction is contact-type V{x) ex 6{x). The critical cur- 
rent Jc, above which stationary current- flow state is ab- 
sent, was found to depend on the strength of barrie r'^^'^" . 

Time-independent GP equation is given by 



2m dx"^ 



*(a;) + [U{x) - pl] *(x) 
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dyV{x~y)\^{y)\H{x)=G. 



(5) 



Substituting a time- independent function ^>{x) = 
yl(a;)e*''''^^) into eq. ([5]), we obtain two equations corre- 
sponding to the real and the imaginary parts 



g[A{x),^{x)]^-^\^A{x)- 
2m dx 



d(f{x) 



n 2 



dx 
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[u{x) - ^A A{x) 



L/2 

dyV{x - y)A{yfA{x) = 0,(6) 
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where A{x) is the amplitude of the condensate wave func- 
tion and (p{x) is the phase of the condensate wave func- 
tion. From eq. ([7]), which is Q for stationary state , it 
follows that 



j^^Aixr^ 

m dx 



const. 
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Substituting eq. ([5]) into eq. (jB]), we can eliminate f{x) 
and obtain the GP equation for the amplitude 

h^ d^ m J2 



^^^^'^^^[-2^d^ + jAW + ^^'^^-^ 

.L/2 ^ 

+ / dyV{x - y)A{y)^ \ A{x) = 0. (9) 

J-L/2 J 

In this paper, we take three kinds of boundary condi- 
tions, as explained below. 

In llllAl and llllBl we use solutions to eq. (O under the 
boundary condition 



A{x + L) = A{x). 

Once we obtain the solution to © with ([1 
sion for (p{x) follows: 

1 
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(10) 
the expres- 
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Here we fix overall phase factor such that 1^9(2; = 0) = 0. 
The values of (p{ztL/2) at the boundaries, which are de- 
termined by A{x) and J through (|TT|) . This set-up de- 
scribes the system attaching a current generator (one can 
regard this one-dimensional system as an orifice attach- 
ing two reservoirs of superfluids or supersolids; See e.g. 
sec. Ill of Ref.E). On the basis of solutions to ([9]) under 
the boundary condition (jlOp . we present the results on 
the phase diagram in IIII Al and Josephson effect in IIIIBI 
In IIII C[ we discuss persistent current in a ring with 
perimeter L. In this case, the boundary condition 



1'{x + L) = *(x) 



(12) 



for (O is used. 

In IIIIDi we discuss nonclassical rotational inertia in 
Hess-Fairbank effect in a ring rotating with angular ve- 
locity uj — 2ttv/L, with use of the solutions to ([5]) satis- 
fying the boundary condition 

^{x + L) = *(x)e-""''^/'"'. (13) 



We can obtain the solutions to eq. ([5]) under p2)) or p^ 
from the solutions to © satisfying ([TU)) by restricting the 
values of J. For the moment, we thus consider the latter 
solutions . 

For consideration of stability of solutions to ([9]) with 
(fTO|) . it is helpful to note, as shown in appendix \K\ that 
^ and ^ result from the stationary conditions of the 
following functional 
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under the constraint ^. Thus the functional K plays the 
role of "energy" under the fixed value of J. The second 
term in eq. (J14p is inherent to the system with constant 
flow as remarked in ref. |38J. 

In the following, we set the potential barrier U{x) to 
be 



Uix) = Uo9id/2~\x\), 



(16) 



where 6{x) denotes the Heaviside step function. V{x^y) 
denotes the soft core two-body interaction:^^'"^'^ 



Vix^y) = Voeia-\x-y\), 



(17) 



where a is the interaction range. 

Let us measure the length, time, energy, ^(x), and 
current in units of a, tq = ma^/h, eo = ?i /ma?', 
y/no = y/N/L, Jq = noh/ma respectively. We define 
a dimensionless parameter g)^ 

2noma^Vo 



9 



(18) 



which characterizes the strength of interaction. 

We solve eq. © in the following procedures: First we 
consider a relaxation dynamics 



dA{x,T) 
dr 



-C[A{x,t)]. 



(19) 



The stationary solutions of eq. p9|) are the solutions of 
eq. (HI). We discretize a; and r with the backward- Euler 
method and then obtain a constituent non-linear equa- 
tion for stationary solutions of eq. ([T^. The resultant 
equation is solved by the Newton method. We perform 
this calculation under a value of /i and repeat it until /i 
satisfies ^. When we obtain different solutions starting 
with different initial conditions for the Newton method, 
we choose the solution with the lower value of energy 
functional K as the lowest energy state for given J. We 
take the system size as L/a = 40, 50. In the following, we 
will present the results for L/a = 50, which are almost 
similar to those for L/a — 40. The number of the mesh 
points in real space is taken as Nj^ = 2^^, 2^^ or 2^^ . We 
take the range of the barrier d to be d/a = 0.5, 1, 1.4, 2. 
We report the results for d/a = 1 in the next section and 
relegated those for d/a = 0.5, 1.4, and 2 to Appendix [Bl 
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FIG. 1. (Color online) Phase diagram in the case of U{x) — 
in the plane of J and g. SF (SS) denotes the superfluid (super- 
solid) phase. In the region with J above the critical current 
(open square), no steady state exists(NS). The solid curve 
represents the Landau critical current in the liquid phase. 
Open triangles represent the liquid-solid phase boundary de- 
termined from the Bragg peak. 
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FIG. 2. (Color online) Density profile for Uo/eo = 0.01 and 
g = 20(< gc), where Jc/Jo = 0.668. 



III. RESULTS 



Phase diagram in the absence of the potential 
barrier 



In the case of U{x) — and J = 0, the property of the 
ground state was studied in Ref. 33; the ground state is 
solid for g > gc = 21.05 • • • and liquid for g < gc and the 
liquid-solid transition at g = gc is continuous. 

Figure [1] shows the phase diagram in the case of 
U{x) — and J > 0. The vertical (horizontal) axis 
is the current J (the interaction strength g). There are 
three regimes: liquid (superfluid) phase denoted by (SF), 
supersolid phase (SS), and non-stationary state (NS). 
The definitions of these regions are as follows. SF is 
the phase in which the lowest energy state for given J 
is |\l'(2;)|/no = 1. SS is the phase in which the lowest 
energy state for given J exhibits a periodic density mod- 
ulation. NS is the state in which no steady flow solution 



of eq. ^ exists for given J. Thus, we can not define the 
phase in this region. Open squares represent the critical 
current. The hquid phase is the region where g < gc and 
J < Ji,. Here Jl = Jl(<?) denotes the "Landau critical 
current" , which is defined by^ 



Jl — no min 



efc 



\hk 



(20) 



in terms of the excitation energy efc in the liquid phase 
with J = 0. In the present model, e^ is obtained analyt- 
ically as^ 



efc 



2m 



2m 



inoVo 



sin ka 



(21) 



with which the boundary (J — Jl) of the liquid phase 
is determined as shown by the blue curve. Open tri- 
angles represent the liquid-solid phase boundary deter- 
mined from the Bragg peak. Those triangles are on the 
blue curve within the numerical accuracy. In Fig. [I] we 
see that for g S [17.5, f^c], the phase changes from the 
liquid phase to the solid phase with J increasing. This 
transition corresponds to the current-driven transition 
between the li quid phase and the modulation phase dis- 
cussed in Refs. |4J, S; andl46l. In the modulation phase, 
the density has the modulation with the period equal to 
the inverse of the roton momentum!^. This phase has 
been discussed to be non-superfluid in Ref. UM and su- 
perfluid in Refs. i45 and 46,. In one-dimensional systems, 
the modulation phase and solid phase belong to a same 
phase, contrary to the systems in higher dimensions. Fig- 
ure [D illustrates how the density profile for g = 20(< ^c) 
changes with J increasing. We confirm that the period of 
modulation coincides with the inverse of roton momen- 
tum and the Bragg peak starts to grow from zero at the 
phase boundary. 

In the excitation energy (|211) . we find that the roton 
minimum exists only for g G [gi{= 9.46 • • •),5c] as noted 
in Ref. [2^. We infer that the Landau critical current 
yields the liquid-solid phase boundary for g > g^ while 
it gives the boundary between the liquid phase and non- 
stationary state for g < g^- Actually, the curves repre- 
senting the critical current and the Landau critical cur- 
rent collapse for g £ [g,:,17.5] as shown in Fig. [T] We 
cannot confirm our expectation numerically. 



B. Josephson effect 

Next we consider the case with a potential barrier. By 
numerical calculation, we find that a steady flow solution 
exists for positive g in the presence of potential barrier. 
Figure [3] shows our numerical results on the spatial de- 
pendences of the density of the condensate (upper panel) 
and the phase of the condensate wavefunction (lower one) 
in the solid phase {g > gc)- As seen in the upper panel, 
the amplitude of modulation becomes larger when g in- 
creases in the solid phase. In the lower panel, the gra- 
dients of the curves correspond to the local velocities of 
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FIG. 3. (Color online) (Upper panel) density profile for 
{g,J/Jo) = (22,0.0921), (25,0.0782), (30,0.0598). For all 
curves, we set Uo/^o = 10. (Lower panel) the profile of the 
phase of condensate wave function for {g, J/ Jo) = (25, 0.0782) 
and t/o/eo = 10~^,C/o/eo = 10. Aip denotes the phase shift 
induced by the barrier. 



the condensate. Small wiggle of ip{x) has the same period 
as that of A{x)'^. This is a consequence of the spatially 
uniformity of the current (j5]). When we compare the 
two curves in the lower panel, we see that the presence 
of the barrier affects the phase of the condensate wave 
function only through the phase shift Aip near x — 0. 
Obviously this is analogous to the DC Josephson effect 
in superconductors^. 

For convenience, the phase-shift Aip is quantified as 



mJ 



L/2 



dx 



L/2 



1 



1 



A{x)^ Ao(x)2 



Aaix) = lim A{x), 



(22) 
(23) 



which is a generalization of the corresponding expressions 
in Refs. Ho and |43 for the liquid phase with superflow. 
Figure |3] shows the relation between J and Aip. In the 
low barrier limit, the J- Aip curve becomes steep for small 
Aip. In the high barrier limit, J- Aip curves approach the 
form 

J = JcSin(A(p), (24) 

which has the same form as the Josephson relation in 




is rewritten as (1101) and 



FIG. 4. (Color online) J-Aip characteristic (the Josephson 
relation) for various value of Uo and for g = 25 > Qc- For a 
technical reason, we take A{x) for C/o/eo = 10~^ as Ao{x) in 
(1221). 




FIG. 5. (Color online) ^/-dependence of the critical current 
Jc for various strength of potential barrier. The solid curve 
represents the Landau critical current. 



superconductorsSA. Comparing Fig.|3]with the Josephson 
relation in Fig. 2 of Ref. |40, one can find that the results 
are very similar with each other for Aip € [0, Aipd with 
Aipc defined by the phase shift at the critical current. 
Figure [5] shows g-dependence of the critical current for 
the Josephson effect. We see that the current-driven solid 
state at g < gc can yield the Josephson effect. In this 
sense, the superflow property of current-driven solid state 
remains even in the presence of an obstacle. 

In Fig. [5l we notice that the J — Jc(.g) curve is not 
affected by the presence of the curve representing the 
Landau critical current for g S [17.5,5c]- We thus expect 
that superflow property does not change so much around 
this boundary for g G [17.5,.gc]- 



C. persistent current in a static ring 

We consider the case where the boundary condition 
(IT^ is imposed. This case describes the condensate form- 
ing a ring with perimeter L. The boundary condition (J12p 



(p{x + L) = if{x) + 27rM, 



(25) 



with an integer M (a winding number). Among the sta- 
tionary solutions obtained in IIIIBl realized are only the 
solutions specified by the discrete set of J satisfying 



L/2 

^ / dy- 



mJ 

'h 



1 



= 2t:M, 



(26) 



-L/2 Mvr 

which immediately follows from (|lip and (|25l) . The so- 
lution specified by a nonzero M, which carries a nonzero 
J, is metastable because the winding number M cannot 
change continuously^; the only way for a current carry- 
ing state to relax to another state specified by a winding 
number M'{^ M) is to cause a phase slip^. When a 
phase slip occurs, the system undergoes an intermediate 
state with a zero point xq of condensate wave function 
^(x) such that ^(xo) — and the phase (p{xo) cannot 
be defined. Such an intermediate state has energy higher 
than the initial state by the local depletion of conden- 
sation energy in the vicinity of a; = xo. The metastable 
state with nonzero M (i.e. nonzero J) yields a persistent 
current if the energy barrier that the metastable state 
has to overcome to relax to another state is sufficiently 
large and the relaxation time is much larger than other 
time scales. 



D. Hess-Fairbank effect, nonclassical rotational 
inertia and superfluid fraction in a moving ring 

On the basis of results in the previous section, we dis- 
cuss supersolidity in the presence of an obstacle moving 
with a constant velocity v under the periodic boundary 
condition (|12p . This situation describes the condensate 
in a toroidal geometry with radius R = L/(27r) and suffi- 
ciently small cross section. The equation under a moving 
obstacle 
,a*(a:,i) ^2 92 , , , 



L/2 



'-L/2 

is rewritten as 

.9^(0;', i) 



dyV{x-y)\'i{y,t)\''^{x,t), 



(27) 



ih- 



dt 



2m dx'^ 

fL/2 



+ [Uix')-fA^'ix\t) 



dy'Vix'^y')\^'iy',t)\H'ix',t),i28) 

L/2 

in terms oi x' = x — vt and ^'(a;',i) — '^{x,t). In the 
moving obstacle, the stable state (we call this state the 
ground state in the moving frame) yields the minimum 
value of the functional^i^ 



£'[^', 4''*] = £■[*', *'*] - vP[ii', *'*], 



(29) 



under the constraint ([2]) and the periodic boundary con- 
dition dni). In dMl), P[^', *'*] denotes the total momen- 
tum 



' dx' 
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2i 






dx' 



dx' 



(30) 



which is nothing but m J_',„ J(x,t)dx with @. The 
ground state in the moving frame is a stationary solution 
of (pS)) because (|^5|) is written as 

rL/2 



i?i 



9*'(a;',i) _ (5£:[*',*' 



<5n;2i*(^'.i)N^' 



(9i 



(5^' 



(5*' 



(31) 



Let the stationary solution to ([T]) for a given J e [— Jc, JJ 
obtained in the previous section be denoted by ^j and 
the accompanying chemical potential fj,j. The stationary 
solution to (|3T]) is then given by 



^'(x') = *j(a;')e™'"^'/'', 



(32) 



with fi = ^j — mv'^/2. Note that from the boundary 
condition (fT2l). it follows that 



(33) 



Substituting ([5^ into ([2^ . we obtain 



(34) 



-E[^j,^j] is even with respect to J and increases with 
J increasing. Thus what we seek is J with minimum \J\ 
such that the boundary condition (155)) is satisfied. From 



ipiL/2) - ^i-L/2) 
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with Aj{y) 
obtain 



L/2 Aj{yy 
^j(y)| follows. From this and (pS)) . we 
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(36) 



with an integer Af . Among the solutions J to ((36|) under 
a given w, let denote by J{v) the solution with minimum 
\J\. Using the results on Aj, we can obtain J{v) as shown 
in Fig. [5] J{v) is periodic with respect to v because v and 
v + 2TTh/{mL) yield the same boundary condition p3p on 
^j. Using this result, we can derive the w-dependence of 
the total momentum, 
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"21 J -L/2 
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FIG. 6. (Color online) w-dependence of J{v) for g — 25 and 
d — a, f/o = eo- Here, vo = fi/mL. 
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L/2 
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* (x ) — tV^ ~ * (^ ) 



dx' 



V / dx\^j{x')\'^ + JL 
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dx' 
Nm,v + mJL 



(37) 



With input shown in Fig. [Bl we plot P as a function of v 
as shown in Fig. [T) The dotted line shows the line P = 
Nmv, which represent the normal state (i.e. the state 
does not exhibit superflow). When \v\ < vq = fi/{mL), 
the slope of the P/{Nrav) ~ v/vq curve represents 1 — 
/s = /n- Here /g = {Nmv — P)/Nmv and /„ represent, 
respectively, the superfluid and non-superfluid fractions. 
Note that Figure [7] can be regarded as a normalized plot 
of the relation between angular momentum L^ — RP and 
angular velocity uj = v / R of the rotating container (ring) 
with the radius R = L/(27r). Compare this figure with 
e.g. Fig. 3.2 of Ref. gg Here wo = vo/R = h/{mLR) 
and /ci — NmR^ (classical rotational inertia). In this 
regard. Figure [7] shows that the supersolid state exhibits 
nonclassical rotation inertia / = /ci/s, which is a sign 
of the Hess-Fairbank effect in supersolid in the circular 
asymmetric container. 



IV. DISCUSSION 

In this paper, we studied one-dimensional supersolids 
in the framework of the mean-field theory. In low- 
dimensional systems, quantum fluctuations are large and 
destroy long-range order^^'--. Therefore, mean-field the- 
ories is not necessarily applicable in low dimensional sys- 
tems. Roughly speaking, our mean-field description is 
valid to describe phenomena over the time scale shorter 
than the inverse of nucleation rates of topological de- 
fects (kink); Those defects destroy the phase coherence 
or positional order in an ordered state. Importance of 
the results in the present paper lies in simplicity of mod- 
els showing the Josephson effect in supersolid phase and 
implication to systems in higher dimensions; Our results 
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FIG. 7. (Color online) u-dependence of total momentum P 
for g = 25 and d — a, Uo ~ eo- Alternatively, this curve can 
be regarded as the relation between angular momentum L^ = 
RP of the condensate and angular velocity uj = v/R of the 
container with the annular geometry with radius R = L/{2tt). 
The unit of angular velocity is luq — h/{mLR). 



strongly suggest that the Josephson effect is likely to oc- 
cur in supersolids in higher dimensions. 

In an earlier work2^, on the other hand, it was con- 
cluded that non-dissipative mass flow around an obsta- 
cle is impossible in a two-dimensional supersolid. There 
are two possibilities of the cause of the difference. One 
is the way of setting boundary conditions. In Ref. |28| . 
the details of the boundary conditions are not avail- 
able. The boundary conditions may determine whether 
or not dissipationless flow is possible. The other possi- 
bility is due to the difference of shapes of the obstacle. 
Although two-dimensional supersolid could exhibit the 
Josephson effect, dissipationless flow circumventing the 
obstacle may be impossible. A two-dimensional solid is 
characterized by positional order and orientational order( 
see e.g. Refs.l53landl54[). A plate-shaped obstacle corre- 
sponding to the Josephson effect may distort only posi- 
tional order, whereas a disk-shaped obstacle may distort 
positional order and orientational order. A future prob- 
lem is to understand these two results, which seemingly 
contradicts with each other, in a consistent way. 

We now discuss how the phase diagram shown in Fig.[T] 
will change in other spatial dimensions where the two- 
body interaction V{r) ex 9{a — r) with a range a{> 0) 
is the same type as that in Eq. ([5]). As in the one- 
dimensional model, the excitation energy e^ in two or 
three dimensions has the roton minimum when g is larger 
than a certain value g^ and the Landau critical current Jl 
becomes zero at a coupling constant 5c(> 5r)- In two or 
three dimensions, the flrst-order phase transition between 
the liquid phase and solid phase at J = was considered 
to occur at a coupling constant g — gi & (5r,ffc)"^^- At 
sufficiently small J, the liquid-solid transition is expected 
to be first-order. The phase boundary between the liq- 
uid phase and a solid phase is represented by a curve 
that starts from {J,g) = {0,gi) and terminates when 
it meets the curve {J,g) = {Jhig),9) representing the 
Landau critical current. Although the modulation phase 



(i.e., the state with laminar pattern of density) in two 
and three dimensions might be no longer stable, there 
still remains a possibility that this phase is metastable 
and contributes as an intermediate state to a transient 
dynamics relaxing from a liquid state to a solid state with 
J>0. 

We also discuss how the critical current Jc{g) of the 
Josephson effect looks like in two and three dimensions. 
We expect the overall feature of the curve J — Jc{g) in 
those dimensions to be similar to that shown in Fig.[5j As 
a notable difference between one dimension and higher di- 
mensions, we expect that the curve J = Jc{g) has a cusp 
or kink when it crosses the first-order phase boundary 
between the liquid and the solid phases. 

From Fig. [3] (lower panel) and Fig. 01 the superflow 
property of the solid phase seems to be similar with 
that of the conventional superfiuid. We discuss a pos- 
sible difference between the two phases. In conven- 
tional Josephson junctions, J-Aip characteristics has a 
curve connecting the points (J, At/?) ~ {Jc,Aipc) and 
{J,Aip) = (0,7r). This curve corresponds to the un- 
stable steady state. In Ref. ^38), it was shown for GP 
equation with contact-type two-body interaction and the 
potential barrier U{x) ex S{x), there exists only a single 
branch that corresponds to the unstable steady state. In 
the present case we study, on the other hand, the num- 
ber of unstable and metastable branches are unknown. 
The structure (i.e., the number and stabilities) of the 
branches corresponding to steady states in the solid phase 
might be different from that in the conventional super- 
fiuid. From the context of non-linear physics, this implies 
that the critical current state in the solid phase in the 
present model might be categorized as a different type 
of bifurcation points from the critical current state in 
the conventional superfiuids. The bifurcation theory dis- 
cusses parameter-dependence of the number and stability 
of stationary solutions of time-dependent non-linear dif- 
ferential equations^^. The bifurcation point means the 
critical parameter at which the number or stability of 
stationary solutions changes. In the present case, equa- 
tion (IT^ has at least one stationary stable solution for 
J < Jc{g) and no stationary solutions for J > Jc(.g). 
The critical parameter J = Jc{g) is thus a bifurcation 
point. It is important to identify the type of bifurca- 
tion points in the present model in the sense that the 
type of the bifurcation in the critical current state gov- 
erns the dynamics in the breakdown of superfiuids, such 
as the emission rate of the topological defects (solitons 
and vortices) around obstacles at J slightly larger than 
the critical valued"—. If the critical current state in the 
solid phase in the present model is identified with a dif- 
ferent type of bifurcation point from that in conventional 
superfiuids, then the dynamical properties in breakdown 
process of superfiuid properties is also expected to be un- 
conventional. The dynamics of the solid phase above the 
critical current and identification of the bifurcation type 
are addressed in a future study. 



V. CONCLUSION 

In conclusion, we showed that superflow properties 
(Josephson effect, Persistent current and Hess-Fairbank 
effect) of a sohd phase are maintained even in the pres- 
ence of an obstacle in a one-dimensional GP model with 
finite-range interaction. A generalization to higher di- 
mensions and the study of the mechanism of breakdown 
of supersolidity are important future issues. 
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Appendix A 

In this appendix, we show that the stationary condition 
of K with respect to the variation A{x) — ?> A{x) + 6A{x) 
and ip{x) — ?> 'fiix) + S(p{x) under the constraint 1^ and 
the boundary condition ^TU\\ . Following the standard pro- 
cedure of variation under a constraint, we consider the 
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L/2 



dx\^{x) 



variation of the functional K ~ K ~ ^ J 
For the functional 

.L/2 

E = E-fi dx\^{x)\^, 

J-L/2 

we obtain 

E[A{x) -^ A{x) + SA{x), ip{x) -^ (p{x) + S(p{x)] 
-E[A{x),ip{x)] 

= SE + OiidA{x)f ,dA{x)Sipix), i5ip{x)f), 
with 



(Al) 



6E = 



dxdA{x)g[A{x),ip{x)] 

-L/2 

dxdip{x)-j-- 

'L/2 



dx 



m 



-Aixf^ 
~" dx 



m 
rn 



dAjx) 
dx 



-1 x=L/2 



SA{x) 



= -L/2 



A{xr^ 

dx 



Sip{x) 



x=L/2 
x=-Lj2 



(A2) 

The first and second terms in RHS vanish when A{x) and 
ip(x) satisfy © and ©. The third term in RHS of (|X2|) 
becomes zero from the boundary condition (|10p . The last 
term becomes 

-1 x=L/2 
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dx 
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FIG. 8. (Color online) d-dependence of the Josephson relation 
for Uo/eo — 10 and g = 25. 
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FIG. 9. (Color online) Density profile for g = 25, J = 0, t/o = 
lOeo, and d — 2a. 



It thus follows that 

SK ^SE- hJ [Stp{L/2) ~ 5lp{-~L/2)] = 0. (A4) 

Appendix B 

In this appendix, we discuss the d-dependence of our 
results {d denotes the range of potential barrier defined in 
(|T6| . Figure [8] shows that the Josephson relation for var- 
ious value of d. We find that d-dependence of the critical 
current is nonmonotonic in contrast to t/o-dependence. 
This behavior of the critical current is due to the com- 
mensurability of the width of the potential barrier and 
the period of density modulations. Fig. IH] shows that the 
density profile for d/a — 2, f/o/eo = 10, J = 0, and 
g ^ 2b. There are two local minima of the density in 
the potential barrier. This result is different from that 
of d < 1.4a. However, Josephson effect occurs for d — 2a 
as shown in Fig. [51 This result suggests that Josephson 
effect in one-dimensional supersolids occurs regardless of 
the width of the obstacle. 
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